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Abstract 

We discuss the relativistic nuclear equation of state (EOS) using a relativistic 
transport model in heavy-ion collisions. From the baryon flow for Au + Au 
systems at SIS to AGS energies and above we find that the strength of the vec- 
tor potential has to be reduced moderately at high density or at high relative 
momenta to describe the flow data at 1-10 A GeV. We use the same dynam- 
ical model to calculate the nuclear EOS and then employ this to calculate 
the gross structure of the neutron star considering the core to be composed of 
neutrons with an admixture of protons, electrons, muons, sigmas and lambdas 
at zero temperature. We then discuss these gross properties of neutron stars 
such as maximum mass and radius in contrast to the observational values. 
PACS number(s): 26.60. +c, 21.65.+f 
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I. INTRODUCTION 



The nuclear equation of state (EOS) at high density is still an unresolved issue though 
many theoretical and experimental efforts have been made in the last two decades to address 
this question in a more systematic way. Theoretically, specially in astrophysics, the den- 
sity of the core inside the compact objects like neutron star is greater than normal nuclear 
matter density, composed of many non-strange and strange degrees of freedom. One of the 
most important characteristic feature of a neutron star is its maximum allowed mass. The 
determination of maximum mass and radius of neutron star are dominated by the interac- 
tions between particles at high density and its EOS. There are many models available in 
the literature to deal with maximum masses of neutron stars. These are relativistic and 
non- relativistic approaches. The non- relativistic models []T],0] based on the potential ap- 
proach describe the nuclear structure for light nuclei. However, the relativistic models 
constructed from Lagrangian approach explain the nuclear structure data for heavy nuclei 
without violating the properties of nuclear matter at saturation density. In both conven- 
tional approaches in the neutron star matter, the estimated maximum masses of neutron 
star are above 2Mq. Recently from several calculations, it has been pointed out p|-p!0| that 
the nuclear EOS should be soft at high density. This is due to fact that all the measured neu- 
tron star masses are less than 2M0 |n|. Various scenarios including the reduced strength 



of vector field, the presence of hyperons and possibility of kaon condensation, have been 
proposed to be soften the EOS. 

Regarding the composition of neutron star matter, there are calculations ||12[ , which in- 
clude kaons as the strange particles along with neutrons and protons e.g., the possibility of 
kaon condensation. Also there are models in the neutron star matter where the compo- 
sition of particles are sigmas and lambdas as strange particles besides neutrons, protons and 
electrons as non- strange particles. Both these proposed models in the neutron star matter 
lead to a soft EOS at high density. In this paper we consider the existence of hyperons in 
the neutron star matter with the recent compiled information of nuclear interaction from 
heavy- ion collisions. 

Experimentally, the nuclear EOS is very important to understand the non- equilibrium 
complicated heavy- ion collisions data at very high energies. Very recently ]13|,|14|, the 
heavy- ion collisions data such as the sideward and elliptic flow have been measured at 
AGS energies. The sideward flow data are mainly determined by the nature of the nuclear 
force in the nuclear EOS. Moreover, the nuclear EOS can be understood better from the 
elliptic flow than sideward flow, because the elliptical flow plays less uncertainties role in the 
opposing stream of matter moving past each other within the reaction plane in the heavy- 
ion collisions. Recently, the beam energy dependence of flow data JT3| , P^ indicate that the 
nuclear EOS is rather soft to lead a possible phase transition to quark gluon plasma at high 
density and hence the strength of the repulsive vector potential must be low to describe 
these data in the heavy- ion collisions. 

In the present discussion, we use an extended version of relativistic mean field model 
including the momenta dependent forces, which are taken into account phenomenologically 
in the relativistic transport model in heavy- ion collisions. We calculate the nuclear EOS by 
using the same dynamic momentum dependence constraints in the nuclear potentials and 
then employ to the neutron star structure calculations. The aim of this paper is to derive 
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the nuclear force from the heavy- ion colhsions data e.g., from nucleon flow data and then 
to study this force on the gross structure of the neutron star by giving less importance to 
the composition of neutron star matter. As far as the strange particles are concerned, we 
take minimum strange particles (E and A) in the neutron star matter calculation at high 
densities. 

The paper is organized as follows: In section 2 we briefly describe the relativistic nuclear 
EOS and its derivation from heavy- ion reactions. In section 3 we employ the same nuclear 
EOS to the neutron star structure with the systematic results. The conclusion and summary 
are presented in section 4. 



II. RELATIVISTIC NUCLEAR EQUATION OF STATE 

The relativistic mean- field theory is very successful model in the relativistic transport 
model in heavy- ion collisions as well as in nuclear structure physics. Originally, Walecka 
had proposed the relativistic mean- field model and later the modified version of this 
has been used widely to calculate nuclear structure and nuclear matter properties. The 
extended version of the Walecka model, so called non- linear relativistic mean- field model 
IQJ^ has the interaction of Dirac nucleons with scalar and vector mesons as well as non- 
linear self- interaction of the scalar field. The extra non- linear self- interaction scalar field 
helps to get the empirical values of bulk properties of nuclear matter at saturation density, 
e.g., the nuclear incompressibility and the value of effective nucleon mass in the desirable 
range. The physics behind this phenomenological successful model is that the nucleon- 
nucleon interaction in the mean- field theory contains strong attractive Lorentz scalar and 
repulsive Lorentz vector components, which almost cancel for low momenta, but produce a 
strong spin- orbit force consistent with the observed single- particle spectra. In the original 
Walecka model the vector potential increases linearly with density, whereas the scalar 
potential changes non- linearly. This is because the vector and scalar potentials have linear 
and non- linear function of density respectively. However, from the heavy- ion collisions data, 
we find that the vector potential also should have non- linear function of baryon density, 
i.e., the strength of the vector potential should be low at high density |]TB|] compared to 
the original model Recently, this fact has been taken by adding the non- linear vector 
meson terms in the original Lagrangian density and applied to the nuclear matter, neutron 
star matter |jTO| and nuclear structure calculations. In our calculation, we take the non- 



linear effect in the vector meson with density by employing the phenomenological momentum 
dependent cut- off to the vector potential term. We adopt this method keeping in mind to 
describe the heavy- ion reactions data at high energies, which generates the nuclear matter 
like situation in the laboratory. We recall that the mean-field energy density for nuclear 
matter in the relativistic mean-field model can be written as Q 

e{m*,nb) = gvVorib - ^m^Kf + - m*Y + -^(m - m*f 

where m* = m — QsSq is the effective nucleon mass, nt is the baryon density and the spin 
and isospin degeneracy is 7 = 4. 5*0 and Vq are the scalar and vector fields with mass 
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and m^, which couple to nucleons with couphng constants gs and g^, respectively. B and 
C are constant parameters describing the scalar self-interactions field and p is the nucleon 
momentum integrated up to the Fermi momentum kf. In (1), the vector and scalar potentials 
depend on density, however, the vector potential increases linearly with density (n^). The 
parameters g^, gs, B and C in (1) are determined by fitting the saturation density, binding 
energy, effective nucleon mass and the compression modulus at nuclear matter density (cf. 
NL3 parameters set from Table I in Ref. 

In our present calculation, we have extended (1) to include a non-linear dependence of 
the vector potential on the baryon density by implementing the momentum (p) dependent 
form factor at the vertices and can be written as 

^O(P) = ^0^^, (2) 

where the cut- off parameters A„i = 0.37 GeV, At,2 = 0.9 GeV and Vq is vector potential. For 
completeness, we incorporate the momentum dependent form factor at the scalar vertices 
in the form given as [in 



P ^^sl 



ys{p) = Vs'-r-i^, (3) 



p2 + A 



s2 



where the cut- off parameters A^i = 0.71 GeV, = 1.0 GeV and Vs is scalar potential. The 
choice of these form factors are similar to that used in effective meson-exchange interactions 
for nucleon-nucleon scattering [|l^ and later this strategy was used in relativistic approach for 



nucleus- nucleus collisions from SIS to SPS energies [0. The values of cut-off parameters 
in vector and scalar vertices are chosen to describe properly the Schrodinger-equivalent 
potential until 1 GeV and the flow data at AGS energies. These cut-off parameters are not 
unique for various type of equation of states to flt Schrodinger-equivalent potential until 1 
GeV and the flow data at AGS energies simultaneously. We note that the form factor eq. (2) 
will make the vector interaction weak at high baryon density and at high energies in heavy- 
ion collisions. At these energies, it has also been observed that the strength of repulsive 
vector potential should be reduced considerably at high density or at high relative momenta 
to describe the flow data. Theoretically, it is important to understand the decrease of vector 
coupling at high density. Contrast to heavy-ion reactions, in this line some works have been 
performed ||19[ and more to be required in details ||20[ . 



We show in Fig. 1, the scalar and vector potential energies as a function of baryon 
density. The solid lines (NLE curves) are associated with the momentum dependent form 
factors given in eqs(2,3) that describe the flow data best from SIS to AGS energies. The 
dashed lines (NL3 curves) are without momentum dependent potentials. The vector part 
for NLE is substantially low at high baryon density as that of NL3 parameter set. At high 
density the reduction of vector potential is more signiflcant than the scalar potential for NLE 
curves. Therefore, the net effect of changed potentials is vector potential due to substantial 
reduction of vector part at high baryon density. For example, at p = Spo, the values of 
vector part and scalar part are 1250 MeV and -511 MeV respectively for NLE, where the 
value of vector part is 1740 MeV and scalar part is -735 MeV for NL3. So, the net reduction 
is dominated by the vector potential in NLE model. The corresponding EOS versus baryon 
density are shown in Fig. 2 for extended momentum dependent model (NLE) as well as the 
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original non- linear model (NL3). NLE has the momentum dependent form factor in the 
vector and scalar potentials. The other nuclear EOS has been discussed more details in the 
Ref. 1^ by varying the nuclear incompressibility from low (soft) 250 MeV to high (stiff) 350 
MeV values. We do not elaborate on that issue here, because we would like to emphasize 
more on the momentum dependent force in the nuclear EOS along the line of heavy-ion 
reaction data. We see in Fig. 2 that NLE nuclear EOS is softer than NL3 at density > 7po 
and is slightly stiffer at density < 7po- The incompressibility is close (~ 380 MeV) to NL3 
value at saturation density. In the next section, we would like to implement this model in 
the neutron star matter, where the core density is in the range of > 5 — 8po- So, in the 
present model, the stiffness of equation of state changes around that density, due to the 
main contribution coming from the reduced vector potential. However, in the heavy- ion 
flow calculation at AGS energies, the stiffness of equation of state not only comes from the 
net reduction of vector potential but also from the transition from hadron to string degrees 
of freedom as discussed in our recent work \T^\. It has been pointed out recently |2T] from 



the simulation calculation that one might even reach lOpo, although only for a very short 
time of a few fm/c at energy range between the AGS and the SPS energies. Hence at AGS 
energy range, the baryon density is expected to reach more than > 5po- 

Recently, the elliptic flow and the sideward flow have been studied theoretically with 
increasing beam energy by various type of equation of state and the possible signature of 
phase transition p2[. More precisely, the beam energy dependence of the observed elliptic 



flow has been interpreted such a possible phase transition. The reason is that a simulation 
model including different kinds of equation of state is consistent with a softening of the 
equation of state. This softening of equation of state can be realized in many ways, for 
example: (i) by reducing the strong repulsive force in the equation of state with help of 
momentum dependent form factor and fitting it with Schrodinger-equivalent potential and 
(ii) by implementing transition from hadronic to string degrees of freedom with beam energies 
in the simulation model []T3[. In our calculation, we implement the former one, where the 



thermodynamic pressure in the extended model NLE is lower as compared to the NL3 model 
due to less repulsive force at AGS energy regime. We thus get reduced repulsive force because 
of the strong cut-off parameters eq. (2) in the vector potential. Also this cut-off makes the 
vector potential non-linear function of baryon density. 



III. NEUTRON STAR MATTER AND PROPERTIES OF NEUTRON STAR 

A. Neutron star matter 

The core of the neutron star plays a significant role to determine the gross structural 
properties like maximum mass and radius of the neutron star. The density of the core inside 
the neutron star is greater than the normal nuclear matter density and hence the nuclear 
interactions are important to construct the neutron star matter EOS around that density. 
Moreover, in such a high density, the strange particles are expected to be present along with 
usual neutron matters like neutrons, protons and electrons. So, in our neutron star matter 
calculation we assume that the core of the neutron star matter is composed of neutrons with 
an admixture of protons, electrons, muons and hyperons (A and S^) [0. The concentrations 
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of each particle can be determined by using the condition of equihbrium under the weak 
interactions (assuming that neutrinos are not degenerate) and the electric charge neutrality: 



f^P l^n /^E! /^A /^n; 
/^S" /^n ~l~ f^ei f^fi /^e; 



P 



ne + n^ + riY.- . (4) 



In addition, the total baryon density is ns = Un + Up + n\ + tiy,- and the baryon chemical 
potential is /i^ = /i„, where rii and /x, stand for number density and chemical potential for 
z— th particle respectively. 

Since the nuclear force is known to favor isospin symmetry and the symmetry energy 
arising solely from the Fermi energy is known to be inadequate to account for the empirical 
value of the symmetry energy (~ 32 MeV), we include the interaction due to isospin triplet 
fh meson in the relativistic non- linear mean field model for the purpose of describing the 
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neutron- rich matter [0. It is noted that the p- meson will contribute a term = g^(np — n„)^ 
to the energy density and pressure. We fix the coupling constant Qp by requiring that the 
symmetric energy coefficient correspond to the empirical value 32 MeV. Then the neutron 
star matter EOS is calculated from energy density e and pressure P are given as follow 0: 

+ y 'S'o + ^ Efg + XI ^FG 

^ i I 

P = \n'X + - l^lSl - ^Sl 

-^St + Y.PFG + T.PFG (5) 

^ i I 

where po is the third component in isospin space. In the above equations Spc and Ppc are 
the relativistic non- interacting energy density and pressure of the baryons {%) and leptons 
(/) respectively. 

The three coupling constant parameters of hyperon- meson interaction are not well 
known. Therefore, we fix the ratio of hyperon-meson and nucleon-meson couplings for cr, 
oj and p mesons respectively (i) by choosing very close to the quark counting rule e.g., 



the potentials seen by A and S in nuclear matter are ~ —30 MeV |^ and (ii) assuming 
the attractive potential seen by A and repulsive potential seen by S to be ~ -30 MeV ||2^ 
and ~ -|-10 MeV P J2^ , |25[| respectively, at nuclear matter density. The analysis of various 
experimental data on hypernuclei p3|-p6| suggest that the strength of S potential may be 
either repulsive or attractive at nuclear matter density. This point will be cleared further 
after analysis of more hypernuclei data in near future and the general discussions are given in 



recent Ref. ||2^. Due to this fact, we consider the two possibility of strength on E potential 



as discussed above. 

Taking all these parameters into the equations (4), we show the concentration of particles 
{xi = rii/uB, i=p, S and A) versus baryon density for NLEl, NLE2 and NL3 models in Fig. 
3. We display p, S and A particles in this figure due to practical importance in neutron 
stars, for example, p fraction plays role for cooling process of neutron stars and the order 
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of appearance of strange particles with density may influence the EOS of the neutron star 
matter. In NLEl and NL3 models, the potentials for A and S are taken to be equal to 
~ -30 MeV, where the potentials for A and S are chosen to be ~ -30 MeV and ~ -|-10 
MeV respectively in NLE2 model. However, the momentum dependent cut-off to the vector 
potential are incorporated in both NLEl and NLE2 models. We notice in Fig. 3 that the 
concentration of particles like S~ and A start appearing after 2 times nuclear matter density 
for all models. In NLEl and NL3 models, the order of appearance of strange particles (first 
E and then A) are same due to equal strength of potential felt by strange particles. Where 
the situation is quite different in case of NLE2 model, here A appears first around > 2.5 
time nuclear matter density and starts coming much later around > 3.5 times nuclear 
matter density. This is due to fact that S sees extra strength +40 MeV potential than 
A potential, which is repulsive. In both NLEl and NLE2 models, the strange particles 
start coming slightly later than NL3, due to reduction of vector potential by momentum 
dependent cut-off as given in equation (2). However, the change of proton concentration 
is not very significant with density for all models, except the slight decreasing tendency at 
high density was shown by NLE2. At around 1.5 times nuclear matter density, the value 
of protons concentration crosses the threshold value 0.11 (horizontal line in Fig. 3), which 
shows that the direct URCA process is possible to lead for cooling of the neutron stars in 



all models p8 



B. Maximum mass and radius of neutron star 

The gross structure of the neutron stars such as mass and radius are calculated from 
the equations that describe the hydrostatic equilibrium of degenerate stars without rotation 
in general relativity is called Tolman-Oppenheimer- Volkoff (TOV) equations [0. From the 
dynamics and transport properties of pulsars, the additional structure parameters of neutron 
stars like the moment of inertia / and the surface red shift z = —, = 1 are important 

and are given more elaborately in Ref. [0. 

We solve the TOV equations by constructing the EOS for the entire density region 
starting from the higher density at the center to the surface density. The composite EOS for 
the entire neutron star density span was constructed by joining the NLE and NL3 neutron 
star matter EOS to that EOS of the density range (i) 10^'* to 5 x 10^° g cm~^ p9|, (ii) 



5 X 10^*^ g cm ^ to 10^ g cm ^ [RH] and (iii) less than 10^ g cm ^ PTI]. The composite 



neutron star matter EOS are plotted in Fig. 4 for NLEl, NLE2 and NL3 models, which are 
used to calculate the neutron star structures as discussed above. From Fig. 4 we find that 
the pressure is low at high density for NLEl, NLE2 EOS and hence the soft EOS compare 
to NL3 EOS. If we look at the Fig. 3, the order of appearance of A particles with density 
reflect in the same order of the nature of EOS. That is, NL3 is stiffer than NLEl and NLE2 
EOS, because the momentum dependent form factor in later two models has reduced the 
vector potential at high density. So, NLE2 is similar to NLEl EOS, except a slight stiffer 
than NLEl due to the strong repulsive potential present on S particles as can be seen in 
Fig. 4. We also notice in Fig. 4 that the NLE2 EOS does not change significantly on the 
choice of repulsive S potential in contrast to NLEl EOS. 

The predicted maximum neutron star masses are very close to the observational values for 



7 



NLEl and NLE2 EOS. The results for the neutron star structure parameters are tabulated in 
Table I and central density vs mass are plotted in Fig 5. From Fig. 5 and Table I, we observe 
that the maximum masses of the stable neutron stars are 2.I8M0, 1.94M0 and 1.97Mq 
and corresponding radii are 11.9km, 10.7km and 10.8km for NL3, NLEl and NLE2 EOS 
respectively. The corresponding central densities are 2.0 x 10^^ g cm~^, 2.2 x 10^^ g cm~^ 
and 2.2 x 10^^ g cm~^ ( > 7 times nuclear matter density) for NL3, NLEl and NLE2 
respectively at the maximum neutron star masses. These maximum masses calculate in our 
models are in the range of recent observations ||3^-p5| , where the observational consequence 



are discussed below. Very recently, it has been observed that the best determined neutron 
star masses |Tl| are found in binary pulsars and are all lie in the range 1.35ibO.O4M0 except 
for the non-relativistic pulsars PSR J1012+5307 of mass M = (2.1 ± O.8)M0 H. There 
are several X-ray binary masses have been measured, the heaviest among them are Vela 
X-1 with M = (1.9 ± O.2)M0 [||] 



34 . From 



and Cygnus X-2 with M = (1.8 ± O.4)M0 
recent discovery of high-frequency brightness oscillations in low-mass X- ray binaries, the 
large masses of the neutron star in QPOAU1820 — 30(M = 2.3)M0 [^] is confirmed and 
this provides a new method to determine the masses and radii of the neutron stars. We also 
tabulate the moment of inertia and the surface red shift in Table I, which are important for 
the dynamic and transport properties of pulsars. 

At this point, we argue that the softening of EOS may lead to kaon condensation in 
neutron stars |^ and hence may give a constraint on the best determined maximum mass 
rT| . However, we feel that from the KaoS data on kaon production, together with kaon 
flow from heavy- ion reactions [0, it is important to know the momentum dependent 
and K~ potentials in dense matter in contrast to the prediction of the chiral perturbation 
theory. In the present calculation, we don't explore this, but work is in progress |38| by 
implementing the same momentum forces as given in eqs (2,3). 



IV. SUMMARY 

We have described the nuclear EOS in the frame work of relativistic mean field theory 
using a relativistic transport model in the heavy-ion collisions. From the heavy-ion collisions 
data, more specifically, the baryon flow for Au + Au systems at SIS to AGS energies and 
above we noticed that the strength of the vector potential has to be reduced substantially 
at high density and high relative momenta to describe the experimentally observed flow 
data at 1-10 A GeV. In a different way, the vector potential should be non- linear function 
of the baryon density. We took this effect by introducing the momentum dependent cut- 
off into the vector potential in contrast to the heavy-ion collision data. We use the same 
dynamic treatment in our relativistic mean field model to calculate the nuclear EOS. It is 
found that the derived nuclear EOS is moderately soft at density > 7po than the original 
considered nuclear EOS without momentum dependent potentials. This due to the reduction 
of repulsive nuclear interaction in the nuclear EOS at high density. We then employ the 
same nuclear EOS to the neutron star structure calculation. 

In the neutron star matter, the core of the neutron stars are considered to be composed 
of neutrons along with an admixture of protons, electrons, muons and hyperons at zero 
temperature. The resulting maximum mass of stable neutron stars are 2.I8M0, 1.94M0 
and I.97M0 for the NL3, NLEl and NLE2 models, respectively. We observed that the 
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maximum mass of the neutron star for NLEl and NLE2 are lower than that for NL3 due to 
a reduction of the vector field at higher densities. Also, we noticed that the potential felt 
by S particles is not so relevant to neutron star structure calculation. The corresponding 
neutron star radii are 11.9 km, 10.7km and 10.8A;m for NL3, NLEl and NLE2, respectively, 
whereas the corresponding central densities are 2.0 x 10^^ g cm^"^, 2.2 x 10^^ g cm~^ and 
2.2 X 10^^ g cm~^ respectively for NL3, NLEl and NLE2 at the maximum neutron star mass. 
We found that the maximum mass for NLEl and NLE2 are in the observable region P^-|55 
1.4Mq < Mffiax < 2.2Mq and the corresponding radius is in between 8-12 km. 



The author likes to thank W. Gassing and A. Ohnishi for critical reading and W. Gassing, 
A. Ohnishi and Y. Akaishi for fruitful discussions. PKS likes to acknowledge the support 
from the Japan Society for the Promotion of Science, Japan. This work is dedicated to Prof. 
Bhaskar Datta, who passed away on 3rd November 1999. 
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Fig.l potential energy per nucleon vs baryon density in units of uq. The solid 
lines (NLE) are momentum dependent potentials and the dashed lines (NL3) 
are without momentum dependent potentials (see the text). 
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Fig. 2 Energy per nucleon vs baryon density in units of no. The models are same 
as Fig. 1. 
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Fig. 3 The concentration of each particle {xi = ni/ns) vs baryon density in units 
of Tif). The momentum dependent potentials have been incorporated in NLEl 
(solid line) and NLE2 (dashed line). The dashed-dot lines (NL3) £tre without 
momentum dependent potentials. The potential are seen by A and E cire same 
in NLEl, NL3 (dashed-dot line) and are different in NLE2 (see text). 
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Fig. 4 The neutron star matter pressure vs the energy density. The models are 
same as Fig. 3. 
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Fig. 5 The neutron star mass vs radius. The models are same as Fig. 3 
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